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ABSTRACT 


Hesperocyparis (Cupressus) goveniana and putative 
Cupressus goveniana subsp. gibsonensis plus H. pygmaea were 
analyzed by Inter-Simple Sequence Repeats (ISSRs). The ISSRs 
analyses support the continued recognition H. pygmaea and H. 
goveniana, but the recognition of Cupressus goveniana subsp. 
gibsonensis Silba was not supported. Phytologia 91(2):277-286 
(August, 2009). 
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The coastal California species of Hesperocyparis (= 
Cupressus) are often found in very small populations and several are 
endangered. The taxonomy of the western hemisphere cypresses has 
recently changed as DNA sequencing of the classical Cupressus 
species (Little et al., 2004; Little, 2006) demonstrated that the western 
hemisphere cypresses are a separate genus from the eastern hemisphere 
cypresses. The eastern hemisphere cypresses maintained the name 
Cupressus and Callitropsis was applied to the western hemisphere 
cypresses plus Callitropsis nootkatensis (D. Don) Oersted ex D. P. 
Little. However, Debreczy, et al. (2009) argued, on morphological 
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grounds, that Callitropsis nootkatensis is a monotypic genus. 
Sequencing of two additional nuclear genes and petN-psbM (Adams et 
al., 2009) supported the recognition of Callitropsis nootkatensis as a 
monotypic genus. Thus, Callitropsis could not be applied to the 
western hemisphere cypresses, so Adams et al. (2009) erected a new 
genus, Hesperocyparis, for these cypresses. 


Silba (2003) proposed the recognition of several new 
subspecies of Cupressus abramsiana Wolf: C. a. subsp. locatellii 
Silba, Eagle Rock, CA; C. a. subsp. opleri Silba, Bracken Brae Forest, 
Santa Cruz, CA; C. a. subsp. neolomondensis Silba, Majors Creek, CA; 
and C. a. subsp. bittanoensis Silba, Butano Ridge, CA. Cupressus a. 
subsp. butanoensis was recognized (Adams et al., 2009) as 
Hesperocyparis abramsiana var. butanoensis (Silba) Bartel & R. P. 
Adams. Silba (2003) also proposed C. goveniana Gordon subsp. 
gibsonensis Silba and C. macrocarpa Hartw. subsp. lobosensis Silba. 
However, the proposed new subspecies are morphologically rather 
indistinct. If these new subspecies are accepted, then these taxa need 
to be considered under the endangered species act. 


Adams and Bartel (2009) examined the volatile leaf oils of H. 
goveniana (Gordon) Bartel, C. g. subsp. gibsonensis Silba and H. 
pygmaea (Lemmon) Bartel (Fig. 1). The leaf oils of these taxa appear 
to separate them into 2 groups (Fig. 1) composed of H. pygmaea and 
goveniana - gibsonensis. There was some variation among the three 
populations of H. pygmaea (Fig. 1) and a hint of differences between 
goveniana and gibsonensis. These differences appear to be in the 
nature of geographic variation. Adams and Bartel (2009) concluded 
that, based on the leaf terpenoids, there was insufficient evidence to 
support the recognition of Silba's Cupressus goveniana subsp. 
gibsonensis. 


To gather additional genetic information about the validity of 
these subspecies, analyses using Inter-Simple Sequence Repeats 
(ISSRs) were conducted. It should be noted that the leaf samples 
utilized in the present study were taken from the same trees analyzed by 
Adams and Bartel (2009). 
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Figure 1. PCO of Hesperocyparis goveniana, Cupressus g. subsp. 
gibsonensis, H. pygmaea plus three plants of C. abramsiana subsp. 
neolomondensis (from Adams and Bartel, 2009). 


MATERIALS AND METHODS 


Plant material - Specimens used in this study: H. abramsiana 
(C. B. Wolf) Bartel, Bonny Doon, Santa Cruz Co., CA, Bartel 1598a-e; 
H. abramsiana var. butanoensis, Pescadero Creek County Park, Butano 
Ridge Fire Rd., San Mateo Co., CA, Bartel 1605a-e.; C. abramsiana 
subsp. /ocatellii Silba, Eagle Rock, Santa Cruz Co., CA, Bartel 1599a- 
e; C. abramsiana subsp. neolomondensis Silba, Wilder Ranch State 
Park, Santa Cruz Co., CA, Bartel 1604a-e; C. a. subsp. opleri Silba, 
Bracken Brae, Santa Cruz Co., CA, Bartel 1600a-e; H. goveniana, SFB 
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Botanical Reserve, Monterey Co., CA, Bartel 1596a-e; C. goveniana 
subsp. gibsonensis Silba, Point Lobos Ranch, Monterey Co., CA, Bartel 
1595a-e; H. pygmaea, Albion Ridge, Mendocino Co., CA, Bartel 
160la-e; Little River Airport, Bartel 1602a-e; Casper Little Lake Rd., 
CA, Bartel 1603a-e, C. macrocarpa subsp. lobosensis Silba, Point 
Lobos State Reserve, Allan Memorial Grove, Monterey Co., CA, Bartel 
1593a-e, East Grove, Bartel 1594a-e; H. macrocarpa (Hartw.) Bartel, 
Crocker Grove, 100 mn of 17 Mile Drive and Madre Lane intersection, 
Monterey Co., CA, Bartel 1597a-e. Bartel specimens are help in his 
personal herbarium in Carlsbad, CA. 


One gram (fresh weight) of the foliage was placed in 20 g of 
activated silica gel and transported to the lab, thence stored at -20° C 
until the DNA was extracted. DNA was extracted using the Qiagen 
DNeasy mini kit (Qiagen Inc., Valencia, CA). ISSR primers were 
purchased from the University of British Colombia (5'-3 seq., annealing 
temperature used'); 807: AGA GAG AGA GAG AGA GT (50°C), 808: 
AGA GAG AGA GAG AGA GC (50°C), 811: GAG AGA GAG AGA 
GAG AC (50°C), 812: GAG AGA GAG AGA GAG AA (50°C), 836: 
AGA GAG AGA GAG AGA GYA (54°C), 840: GAG AGA GAG 
AGA GAG AYT (54°C), 841: GAG AGA GAG AGA GAG AYC 
(54°C), 847: CAC ACA CAC ACA CAC ARC (58°C), 861: AGC AGC 
AGC AGC AGC AGC, (58°C), 881: GG TGG GGT GGG GTG (50°C), 
886: VDV CTC TCT CTC TCT CT (50°C), 887: DVD TCT CTC TCT 
CTC TC (54°C), 895: AGA GTT GGT AGC TCT TGA TC (50°C), 
900: ACT TCC CCA CAG GTT AAC ACA (50°C). 


PCR stock solutions (Taq, primer, and buffer) were made in 
bulk so that all the PCR reaction tubes for a primer were prepared using 
the same bulk stock. This is a critical factor for minimizing variation in 
band intensitics from sample to sample (see Adams, Flournoy and 
Pandey, 1998, for protocols to minimize PCR band variation). PCR 
was performed in a volume of 15 wl containing: 7.5 ul Epi-Centre 2X 
buffer E (containing 0.4 mM of each dNTP, final conc. = 0.2 mM), 0.75 
ul primer (0.6 uM final conc.), 0.75 ul Epi-Centre Fail-Safe Taq (0.75 
unit/rxn.), 6 ul genomic DNA (0.3 ng/rxn.). A negative control PCR 
tube containing all components, but no genomic DNA, was run with 
each primer to check for contamination. DNA amplification was 
performed in an MJ Programmable Thermal Cycler (MJ Research, 
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Inc.). Samples were run in duplicate to insure reproducibility (Adams, 
Flournoy and Pandey, 1998). A temperature profile was obtained for 
each well of the thermocycler to be sure that no variation existed 
among wells in the heating/ cooling block. The thermal eyele used 
was: 94° C (1.5 min) for initial strand separation, then 39 eyeles of 91° 
C (1 min), 50°C (or 54°C or 58°C, see above) (2 min), 72° C (2 min). 
Two additional steps were used: 91° C (1 min), 50°C (or 54°C or 58°C) 
(2 min) and 72° C (5 min) for final extension. The temperature inside a 
PCR tube containing 15 ul buffer was monitored with a temperature 
probe, quantitated and printed for each step for each of the 40 cycles for 
every PCR run (Adams, Flournoy and Pandey, 1998) to insure that each 
eycle met temperature specifications and that each PCR run was exactly 
the same. Amplification products were analyzed by electrophoresis on 
1.5% agarose gels, 70V, 55 min, and detected by staining with ethidium 
bromide. The gels were visualized over UV light and scanned to digital 
images. The digital images were size normalized in reference to 
pGem® DNA size markers before band scoring. Bands were scored as 
present (4 = faint, 5 = bright, 6 = v. bright) and absent (0). Bands that 
were inconsistent in replicate analyses were not scored. 


Associational measures were computed using absolute 
character state differences (Manhattan metric), divided by the 
maximum observed value for that character over all taxa (= Gower 
metric, Gower, 1971; Adams, 1975). Principal coordinate analysis 
(PCO) was performed by factoring the associational matrix using the 
formulation of Gower (1966) and Veldman (1967). It should be noted 
that problems of homology of RAPD DNA bands on agarose gels can 
be significant (Rieseberg, 1996), but these errors can be accounted for 
using multivariate statistical methods (PCO) (see Adams and 
Rieseberg, 1998). A minimum spanning diagram was constructed by 
selecting the nearest neighbor for each taxon from the pair-wise 
similarity matrix, then connecting those nearest neighbors as nodes in a 
network (Adams, et al. 2003). 


RESULTS AND DISCUSSION 
The use of 14 ISSR primers resulted in 98 scoreable bands 


among the taxa. PCO of the association matrix removed five 
significant eigenroots accounting for: 15.2, 14.0, 9.1, 7.1, and 6.2% of 
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the variance among the OTUs. It is noteworthy that there was 
considerable variation among individuals and this is reflected in the 
relatively small amount of variance that was extracted by the first 3 
eigenroots. 


Ordination of the taxa shows that H. goveniana and H. 
pygmaea are very well resolved from H. abramsiana and H. 
macrocarpa (Fig. 2). 


There appears to be considerable variation among the H. 
pygmaea individuals (Fig. 2). In addition, there is also variation among 
the goveniana - gibsonensis individuals (Fig. 2). 
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Figure 2. PCO of H. goveniana, Cupressus goveniana subsp. 
gibsonensis, and H. pygmaea with exemplars of H. abramsiana and H. 
macrocarpa. 
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In order to better visualize the variation among populations of 
H. pygmaea, a new PCO was performed which contained only the H. 
pygmaea individuals plus two exemplar Hf. goveniana samples. The 
PCO using 78 ISSRs resulted in four cigenroots accounting for 16.6, 
12.8, 10.0, and 8.9% of the variance among individuals. Clearly there 
are considerable differences among H. pygmaea individuals (Fig. 3). 
There is a trend for the Casper Little Lake population (plus 2 
individuals from the Albion Ridge population) to cluster together (Fig. 
3). The balance of the Albion Ridge and Little River population plants 
are interspersed (Fig. 3). There was also a trend in the terpenoids to 
subdivide these populations (Fig. 1), but the pattern is slightly different 
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Figure 3. PCO based on 78 ISSRs for H. pygmaea. Two individuals of 
H. goveniana were included in the PCO. 
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than in the ISSR data (Fig. 3). It seems likely that these populational 
differences are normal geographic variation and should not to be 
recognized as separate taxa. 


The variation between H. goveniana and C. g. subsp. 
gibsonensis was further examined using 75 ISSRs for their ten samples 
plus two individuals of H. pygmaea. PCO resulted in four eigenroots 
that accounted for 23.1, 19.5, 11.7, and 10.2% of the variance among 
these OTUs. Ordination reveals (Fig. 4) that H. goveniana and putative 
C. g. subsp. gibsonensis form a large group, with some separation 
between the taxa. However, as these putative subspecies are each in 
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Figure 4. PCO of H. goveniana and C. g. subsp. gibsonensis 
individuals, plus two exemplars of H. pygmaea. 
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distinct (and disjunct) populations, geographic differentiation could 
well explain this clustering. 


In summary, taking both the terpenoid and ISSR data into 
consideration, there appears to be sufficient genetic differentiation to 
support the recognition of H. pygmaea, but there is insufficient 
differentiation to support the recognition of Silba's Cupressus 
goveniana subsp. gibsonensis. 
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